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Abstract 

We provide a comparison of the results of four SUSY mass spectrum calculations 
in mSUGRA: Isajet, SuSpect, SoftSusy, and SPheno. In particular, we focus on the 
high tan P and focus point regions, where the differences in the results are known 
to be large. 
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1 Introduction 

Many SUSY studies rely on computer codes that calculate the mass spectrum of the 
minimal supersymmetric standard model (MSSM), the couplings, branching ratios, etc., 
from given sets of model parameters. For the LHC, for instance, many simulations are 
done for particular benchmark scenarios or by mapping the (mg, "^1/2) parameter plane. 
For such studies it is certainly important wheather a particular decay channel is open 
or not and what branching ratio it has. Also, theoretically or experimentally excluded 
regions depend on the details of the spectrum. Studies for an e^e~ Linear Collider deal, in 
addition, with high precision measurements of (s)particle properties, with the determina- 
tion of the underlying SUSY breaking parameters, their extrapolation to the GUT scale, 
model destinction, etc. Experimental accuracies of the per-cent or even per-mille level are 
expected. It is thus clear that we need theoretical predictions of a precision comparable 
to the experimental accuracy. However, it has been noticed that different programs 
can give quite different results for the same set of input parameters. 

In this article, we compare the mass spectrum calculations of four public codes: 
Isajet 7.63 §, SuSpect 2.005 [§, SoftSusy 1.4 §, and SPheno 1.0 §, in the minimal su- 
pergravity (mSUGRA) framework. We discuss the renormalization group (RG) running 
and the implementation of radiative corrections, concentrating on the parameter regions 
where the largest differences are encountered: large tan/5 and large mo. An overview of 
which corrections are implemented in each of the four programs is given in Table 1. 



2 Large tan (3 

Large tan/5 has always been recognized as a difficult case since it requires a thorough 
treatement of the bottom Yukawa coupling hi,. It is well known that /if, gets large 
tan/5 enhanced corrections from SUSY loops, the dominant contributions coming from 
hg and tx^ exchanges. These generate a ifg^fe coupling, which is forbidden at tree-level, 
£ ^ hb H^bb + Ahb -ff^^fe. This modifies the tree-level relation between the bottom mass 
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Table 1: RGEs and radiative corrections implemented in Isajet, SuSpect, SoftSusy, and 
SPheno. 



and Yukawa coupling, nib = hbVi — ^ mf, = hhVi{l + A;,) with = (Ahb/hb) tan/5. In the 
programs under discussion this is taken into account as 

hb{Mz) = mr^^(Mz)M(M^) , mr'^(Mz) = . (1) 

Here rhf^ is the DR bottom mass in the Standard Model and Arrib = (Amf,) 
contains the SUSY-loop corrections. The complete 1-loop expression for Arrib is given in 
1^.0 Compared to the naive 1-loop expansion Tfij^ssM _ ^sm^^]^ — Amb/mb), eq. (|I|) makes 
a numerical difference of about 10% in hb and about 10-30% in for large tan f3. The re- 
summation of SUSY threshold corrections [|l^ will be discussed elsewhere fl^ . Although 
all four programs now apply eq. (|l]), some numerical differences in hb remain. These are 
partly due to differences in a^: Suspect, SoftSusy and SPheno calculate in the DR 
scheme, Isajet uses the MS value. Another reason is that Isajet uses nib = f^biMsusv) 
for the expression Anib/nib in eq. (|l]), while the other programs use nib{Mz) or the bottom 
pole mass; also the gluino masses differ by about 5%. Moreover, the vacuum expectation 
values Vi^2 are not running in Isajet. 

The bottom Yukawa coupling has its largest effect in the Higgs sector. Figure |I| shows 
the running of ^ for tjiq = 400 GeV, mi/2 = 300 GeV, Aq = 0, fi > 0, and the two 
cases tan/3 = 10 and tan/? = 50. As one can see, there is good agreement for not too 
large tan/3. However, for tan/5 = 50, quite different results are obtained for m'jj , whose 
evolution is driven by hb- 



dmjj^ 3 



dt Svr 



hbXb + ..., Xb = {ml + ml + + Al) 



^Here note that |§, |[ D and partly have different conventions, e.g., for the ordering of the squark 
mass eigenstates and the sign of /it. 
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Note in particular the dotted line which shows the result obtained with Isajet 7.58. In 
this version, the SUSY corrections to /if, were not yet resummed. 




Figure 1: Running of m'jj^^ as a function of the scale Q, for mo = 400 GeV, mi/2 = 
300 GeV, Ao = 0, /i > 0, tan/5 = {10,50}, Mt = 175 GeV. The full (dotted) lines are 
for Isajet 7.63 (7.58), the dashed lines are for SuSpect 2.005, the dash-dotted ones for 
SoftSusyl.4, and the dash-dot-dotted ones for SPhenol.O. 

The differences in directly translate into and thus into the physical Higgs 
boson masses, since 

™i = f + + (3) 

C2/3 Vi V2 

Here fnjj. = m'jj_ — ti/vi, i = 1,2, and ti^2 are the tadpole contributions. The self energies 
of Z and A have been neglected in eq. (^. We note that including only the tadpoles 
from the third generation is in general a good approximation. The remaining 1-loop 
contributions account for a C(l%) correction. 

Figure |^ shows the Higgs boson masses obtained by the four programs as a function 
of tan p. The new Isajet version 7.63 has led to a majer improvement compared to 
the situation discussed in [|, |l5l (the results obtained by Isajet 7.58 are again shown as 
dotted lines in Fig. For and there is now agreement within ~ 10% up to 
tan/? ~ 45. Sources for the remaining differences are pointed out above. Moreover, it 
makes a difference whether one uses running couplings and/or masses for the tadpoles 
ti^2- Here each program has a different approach. For the neutral scalars, however, the 
situation is not so good. Especially for m/jO, a discrepancy of ~ 4 GeV is too large 
compared to the expected experimental accuracy. This discrepancy is mainly due to 
the different radiative corrections taken into account for the {h'^, H^) system. They vary 
between 1- and 2-loop, effective potential and diagramatic calculations, see Table 1. Given 
the expected experimental accuracy for rriho it is clear that the best available calculation 
should be used. 
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Figure 2: Higgs boson masses as a function of tan/3, for tjiq = 400 Gev, mi/2 = 300 GeV, 
^0 = 0, /i > 0, Mt = 175 GeV; full (dotted) lines: Isajet7.63 (7.58), dashed: SuS- 
pect 2.005, dash-dotted: SoftSusy 1.4, dash-dot-dotted: SPhenol.O. 

3 Large mo 

For large mo, the running of m^^^ becomes very steep and very sensitive to the top Yukawa 
coupling ht = rht/v2: 

^-^,h,X, + ..., X, = {ml + ml + ml^+Al). (4) 
As a result, the /i parameter given by 



niH^ — rn^j^ tan^ j3 1 
tan2 (3-1 2 



becomes extremely sensitive to ht. This is visualized in Fig. |^ where we show in (a) 
the running of m?^^ ^ for mo = 1450 GeV, and in (b) /x as a function of mo. The other 
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Figure 3: a) Running of m'jj_^ ^ for mo = 1450 GeV; b) // as a function of mo; for mi/2 = 
300 GeV, Ao = 0, /i > 0, tan/3 = 10, and Mt = 175 GeV; full (dotted) lines: Isajet 7.63 
(7.58), dashed: SuSpect 2.005, dash-dotted: SoftSusy 1.4, dash-dot-dotted: SPhenol.O. 



parameters are mi/2 = 300 GeV, Ao = 0, /i > 0, and tan (3 = 10. The large discrepancy 
in n for mo > 800 GeV lead to completely different chargino/neutralino properties and 
likewise to very different excluded regions in Isajet compared to the other programs. For 
instance, radiative EWSB breaks down in Isajet for mo ^1-5 TeV. In SuSpect, SoftSusy, 
and SPheno, this happens only for mo > 2.5-2.8 TeV.0 

In order to understand the behaviour in Fig. it is useful to write eq. (^ in the form 

/i^ ~Cim^ + C2m?/2-0.5M|. (6) 



Approximate analytical expressions for ci and C2 can be found e.g., in [|1^, ^7j. For Aq = 
and tan/3 = 10 we get []T^ 

Since the Higgs potential is minimized at Msusy = we take rht in eq. (|^ as 

rfit = 7ht{MsusY)- The mo dependence seen in Isajet is reproduced for mt ~ 151 GeV. 
The one of SuSpect, SoftSusy and SPheno is reproduced for rht ~ 155 GeV. Figure | 
shows a contour plot of /i in the {mo, rht) plane. Notice the fast increasing dependence 
on rht for increasing mo. Notice also that for hit ~ 156-157 GeV, fi becomes almost 
independent of mo, which is the actual focus point condition. 

There are some obvious differences in the calculations. For instance, Msusy, the scale 
where the SUSY parameters are frozen out and the Higgs potential is minimized, varies 
by about 100 GeV due to different radiative corrections to the stop masses, c.f. Table 1. 
In the loop corrections to mt, analogous differences occur as discussed above for Amfe/m^. 

After the conference, a sign error was corrected in SPheno. As a consequence, its rcsuhs for large 
niQ now nicely agree with those of SoftSusy and SuSpect (contrary to what was presented in the talk). 
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Figure 4: The parameter [i as given by eq. (P) in the (mo, m*) plane, for mi/2 = 300 GeV, 
Aq = 0, and tan (3 = 10. 



Also the evolution of ht between Mz and Mt and the inclusion of threshold effects are 
delicate points. However, this is not yet sufficient to explain the observed discrepancies. 
More work is needed to clarify the situation. 



4 Conclusions 

For the calculation of the SUSY mass spectrum from GUT scale boundary conditions, 
there are two particular difficult parameter regions where large numerical differences have 
been noticed: large tan/3 and large mo. These regions are very sensitive to the bottom 
and top Yukawa couplings, respectively. 

The inclusion of the SUSY 1-loop corrections to hb has led to a considerable improve- 
ment in the large tan/3 case. In particular, the four programs now agree on niA within 
< 10% for tan/? < 45. Further improvements are of course desirable. 

For large mo, there are still very large numerical discrepancies due to the corrections 
to hf. As a matter of fact, hf is much smaller in Isajet than in the other programs. Some 
differences in the calculation of ht have been pointed out, but these do not satisfingly 



explain the observed discrepancies. Work is in progress to clarify the situation [Q. 



Acknowledgements 

We thank H. Baer, A. Djouadi, and J.-L. Kneur for useful discussions. 



References 



6 



[1] 

[2] 

[3] 
[4] 

[5] 

[6] 
[7] 

[8] 
[9] 



B. C. Allanach, in Proc. of the APS/DPF/DPB Summer Study on the Future of Par- 
ticle Physics (Snowmass 2001 ) ed. R. Davidson and C. Quigg, |arXiv:hep-ph/0110227| . 

N. Ghodbane and H. U. Martyn, in Proc. of the APS/DPF/DPB Summer Study 
on the Future of Particle Physics (Snowmass 2001) ed. R. Davidson and C. Quigg, 
^rXiv:hep-ph/020l233| . 



H. Baer, F. E. Paige, S. D. Protopopescu and X. Tata, |arXiv:hep-ph/0Q01086 
pittp : / /paige . home . cern . ch/ paige/| 

A. Djouadi, J.-L. Kneur and G. Moultaka, 

|http : / /www ■ 1pm . univ-montp2 . f r : 6714/^kneur/ suspect . html 



B. C. Allanach, Comput. Phys. Commun. 143 (2002) 305 ||arXiv:hep-ph/0104145|| 



|http : / /allanach . home . cern . ch/ allanach/ sof tsusy . html 



W. Porod, in preparation. 

D. M. Pierce, J. A. Bagger, K. T. Matchev and R. J. Zhang, Nucl. Phys. B 491 
(1997) 3 [|arXiv:hep-ph/9606211|] . 



M. A. Bisset, Ph.D. thesis, Univ. Hawaii, UMI-95-32579, 1995. 



M. Carena, M. Quiros and C. E. Wagner, Nucl. Phys. B 461 (1996) 407 ||arXiv:hep- 
ph/9508343|| . 



[10] 
[11] 



S. Heinemeyer, W. Hollik and G. Weiglein, Comput. Phys. Commun. 124 (2000) 76 
|^rXiv:hep-ph/98l2320|| . 



G. Degrassi, P. Slavich and F. Zwirner, Nucl. Phys. B 611 (2001) 403 ||arXiv:hep^ 
ph/0105096|| ; A. Brignole, G. Degrassi, P. Slavich and F. Zwirner, Nucl. Phys. B 631 



(2002) 195 ||arXiv:hep-ph / 1 1 2 1 77|| , |arXiv: hep-ph / 206 1 1 



[12] 
[13] 

[14] 

[15] 

[16] 
[17] 



R. Hempfling, Phys. Rev. D 49 (1994) 6168; L. J. Hall, R. Rattazzi and U. Sarid, 
Phys. Rev. D 50 (1994) 7048 ||arXiv:hep-ph/9306309|| . 



M. Carena, D. Garcia, U. Nierste and C. E. Wagner, Nucl. Phys. B 577 (2000) 88 
|^rXiv:hep-ph/99l25l6|| . 



B. Allanach, S. Kraml and W. Porod, in preparation. 

S. Kraml, talk at the 2nd workshop of the extended ECFA/DESY study "Physics and 
Detectors for a 90 to 800 GeV Lmear Collider", 12-15 Apr 2002, St. Malo, France. 

M. Carena, M. Olechowski, S. Pokorski and C. E. Wagner, Nucl. Phys. B 426 (1994) 
269 [ |arXiv:hep-ph/ 940225311 . 



M. Drees and S. P. Martin, |arXiv:hep-ph/9504324 . 



7 



